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DESCRIPTION 

Methods for Detecting and 
Evaluating Heart D isorders 

Cross-Raference to Related Appli cation 

This application is a continuation-in-part of 
copending application S.N. 541,881, filed June 20, 1990 in 
the name of the same inventors and with the same title. 

5 Background of the Invention 

1. Field of the Invention 

This invention relates to heart disorders. More 
specifically, this invention relates to detecting and 
evaluating arrythmia, fibrillation and related disorders 
10 by manipulation of an electrocardiogram signal. 

2. Descripti on of Related Art 

Despite major advances in the diagnosis and treatment 
of ischemic heart disease over the past decade, a 
substantial number of patients each year may suffer sudden 

15 cardiac death as a consequence of ventricular fibrillation 
(VF) . To date, no reliable predictive or preventive 
measures have been developed. By all outward appearances, 
VF is a highly complex, seemingly random phenomenon. So 
are other related heart disorders, including those stages 

20 in heart behavior which typically precede VF (onset of 
VF) . Accordingly, it is difficult for automated devices 
to determine with any reliability that a patient is 
undergoing VF or onset of VF. Moreover, onset of VF may 
also be difficult to determine with any reliability, even 

25 for skilled medical personnel. 

A method of detecting and evaluating heart disorders 
would therefore find wide applicability and utility. 
Patient monitoring devices may summon medical personnel if 
the patient is undergoing VF or onset of VF. Automatic 
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devices which attempt to counter VF, e.g. automatic 
implantable cardiac defibrillators (AICDs) may vary their 
operation based on evaluation of the severity of the 
patient 1 s condition. Methods for reliably evaluating the 
5 risk of VF may also have important utility in monitoring 
patients undergoing surgery or other critical therapy* 

It has been found that some anti-arrhythmic drugs may 
also have a pro-arrhythmic effect in excess 
concentrations. For example, guinidine has been known to 

10 be toxic in this manner. A method of detecting and 
evaluating heart disorders would also have wide 
applicability and utility in determining if a patient has 
been subjected to a toxic (or partially toxic) dosage of 
a drug relating to heart condition. 

15 Chaos theory is a recently developed field relating 

to phenomena which appear to be highly complex and 
seemingly random, but which may be described as the 
deterministic result of relatively simple systems. Chaos 
theory may have potentially wide applications in biologic 

20 and other systems involving ambiguity and uncertainty. 
For example, it has been conjectured that chaos theory may 
be valuable for describing certain natural processes, 
including electroencephalogram (EEG) and electrocardiogram 
(EKG) signals. Techniques for detecting and evaluating 

25 aspects of deterministic chaos are known in the field of 
chaos theory, but have found little application in the 
medical field. 

Accordingly, there is a need for improved methods and 
devices for detecting and evaluating heart disorders, 

30 including ventricular fibrillation (VF) and the onset of 
VF. 

Summary of the Invention 

A first aspect of the invention provides a method for 
detecting a heart disorder, by examination of a phase- 
35 plane plot (PPP) of a patient electrocardiogram (EKG) . A 
normal patient will have a PPP which is relatively smooth; 
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a patient at risk of developing ventricular fibrillation 
(VF) onset will have a PPP which exhibits features of a 
chaotic process, such as multiple bands , "forbidden 
zones", periodicity with period-doubling and phase 
5 locking? a patient exhibiting VF will have a PPP which 
appears noisy and irregular. Differing PPPs may be 
readily recognized, thus detecting patients with heart 
disorders . 

In a preferred embodiment, the PPP's degree of 
10 deterministic chaos may be measured by a processor, such 
as by graphic and numeric analysis. (1) The processor may 
measure a Lyapunov exponent or a fractal dimension of the 
PPP. (2) The processor may determine a Poincare section 
of the PPP and examine that Poincare section for 
15 indicators of deterministic chaos. Also, the processed 
PPP arid Poincare sections may be reviewed by a human 
operator. 

A second aspect of the invention provides a method 
for detecting a heart disorder, by examination of a 

20 frequency-domain transform (suph as an FFT) of a patient 
EKG. A normal patient will have an FFT with a discrete 
spectrum, while a patient exhibiting VF will have an FFT 
with a relatively continuous spectrum and a peak energy at 
a relatively low frequency (e.g., about 5-6 Hz). A 

25 patient exhibiting VF which is difficult to revert with 
shock will have an FFT with a peak energy at a relatively 
high frequency (e.g., about 10 Hz or more). 

In a preferred embodiment, an automatic defibrilla- 
ting device may comprise means for delivering a variable 

30 shock, the size of which is determined at least in part by 
the FFT's peak energy. The def ibrillating device may also 
comprise means for signalling an alarm if the FFT's peak 
energy is at a relatively high frequency. 

A third aspect of the invention provides a method for 

35 detecting drug toxicity, based on particulars of an 
action-potential duration (APD) restitution curve, or an 
action-potential amplitude (APA) curve, which is con- 
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structed for the patient, such as fitting an exponential 
relation to that curve or such as a parameter time 
constant for that curve. The slope of the fitted curve 
will indicate the patient 1 s possibility of predisposition 
5 to arrythmia. Differences in the parameters of the fitted 
curve allow one to distinguish between normal and abnormal 
patients, e.g. those at risk of arrythmia or ischemia. A 
normal patient will have a relatively low parameter time 
constant; a patient who is exhibiting drug toxicity will 
10 have a relatively high parameter time constant. A PPP of 
APD or APA data may also be generated, and the analytical 
techniques described herein may be utilized to interpret 
that PPP, to determine and evaluate drug toxicity. 

^yj T ^f pegcrjptjQfr of the prawjtpqs 
15 Figure 1 shows a patient monitoring system. 

Figure 2 shows a set of sample EKG signals. 
Figure 3 shows a set of corresponding PPPs for the 
sample EKG signals of figure 2. 

Figure 4 shows an example PPP and a corresponding 
20 Poincare section. 

Figure 5 shows an example PPP and a corresponding 
time-lapse Poincare section. 

Figure 6 shows a set of corresponding frequency- 
domain transforms, obtained by performing a fast Fourier 
25 transform (FFT) on the EKG signal. 

Figure 7 shows an improved automatic implantable 
cardiac defibrillator ("AICD") . . 

Figure 8 shows a signal response of an individual 
heart muscle cell to a stimulus, known in the art as 
30 "action potential". 

Description of the Preferred Embodiment 

A first aspect of the invention relates to detection 
and evaluation of heart disorders by examination of a 
phase-plane plot (PPP) of a patient electrocardiogram 
35 (EKG) . 
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Figure 1 shows a patient monitoring system. A 
patient 101 is coupled to an electrocardiogram (EKG) 
device 102, which acquires EKG signals and transmits them 
to a processor 103. The processor 103 may display the EKG 
5 signals on a monitor 104 (as is well-known in the art) , or 
it may process the EKG signals and display any results of 
processing on the monitor 104. 

EKG signals are well-known in the art, as are methods 
of acquiring them. As used herein, an EKG refers to a 

10 surface electrocardiogram, but other forms of electro- 
cardiogram would also work with the methods disclosed 
herein, and are within the scope and spirit of the 
invention. For example, the EKG shown herein may comprise 
a surface EKG, an epicardial EKG, an endocardial EKG, or 

15 another related signal (or set of signals) measured in or 
near the heart. Moreover, the signal which is manipulated 
may be a voltage signal, a current signal, or another 
related electromagnetic values (or set of values) . 

Figure 2 shows a set of sample EKG signals. A first 

20 EKG signal 201 shows a normal patient. A second EKG 
signal 202 shows a patient in transition to VF. A third 
EKG signal 203 shows a patient with VF. 

The processor 103 may construct a phase-plane plot 
(PPP) from the EKG signal. A first type of PPP comprises 

25 a plot of an EKG variable against its first derivative. 
In a preferred embodiment, the EKG variable is voltage, v 
(itself a function of time); its first derivative is dv/dt 
(also a function of time) • 

However, it would be clear to one of ordinary skill 

30 in the art, after perusal of the specification, drawings 
and claims herein, that wide latitude in construction of 
the PPP is possible. The variable chosen for the PPP may 
be any one of a variety of different parameters, including 
EKG voltage, current, or another signal value. The chosen 

35 variable (v) may be plotted against its first time 
derivative (dv/dt), its second time derivative d 2 v/dt 2 , or 
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another time derivative d^v/dt 11 . Or, an Mth derivative may 
be plotted against an Nth derivative. 

Another type of PPP may comprise a plot of an EKG 
variable (or an Nth derivative thereof) against a time 
5 delayed version of itself, (e.g. v(t) versus v(t-$t)). 
This type of PPP is sometimes also called a "return map". 
This type of PPP is less sensitive to EKG signal noise. 

Another type of PPP may comprise a plot of three EKG 
variables (or Nth derivatives thereof) simultaneously 
10 (e.g., v, dv/dt, and d 2 v/dt 2 ) . Such a PPP would be 3- 
dimensional. Where the PPP is 3 -dimensional, it may be 
displayed stereoscopically, or a 2-dimensional plane "cut" 
of the 3-dimensional display may be displayed on a 2- 
dimensional display. It would be clear to one of ordinary 
15 skill in the art r that all of these choices described 
herein, or combinations thereof, would be workable, and 
are within the scope and spirit of the invention. 

Figure 3 shows a set of corresponding PPPs for the 
sample EKG signals of figure 2. A first PPP 301 
20 corresponds to the first EKG signal 201. A second PPP 302 
corresponds to the second EKG signal 202. A third PPP 303 
corresponds to the third EKG signal 203. 

Part of this aspect of the invention is the discovery 
that a normal patient will have a PPP which exhibits the 
25 regularity and smoothness of an EKG signal from that 
normal patient, while a patient undergoing VF will have a 
PPP which exhibits the irregularity and complexity of an 
EKG signal which might be deterministic chaos (e.g., a 
periodicity, banding and "forbidden zones"). Moreover, a 
30 patient in transition from normal into VF (i.e., in VF 
onset) exhibits a PPP which is consistent with an 
assessment that the EKG signal for the patient is in 
transition to deterministic chaos. 

A normal patient has a relatively regular beat-to- 
35 beat EKG signal. As the patient transitions to VF, the 
patient's EKG signal at first shows oscillations between 
pairs of alterant regular beat-to-beat signals. As the 
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transition continues, the patient's EKG signal then shows 
oscillations between greater and greater numbers of 
alterant regular signals (e.g., four possible alternants, 
eight possible alternants, etc.)/ until it is no longer 
5 possible to identify alterant regular signals and the EKG 
signal is irregular and highly complex. At that point, 
the patient is generally said to be exhibiting VF. 

In like manner, the patient's PFP will transition 
from a smooth single-banded display, through a multi- 
10 banded display (showing multiple alternants) and finally 
to an irregular and highly complex display. The display 
change in the PPP is so striking that even a relatively 
untrained person can see the difference. This is in 
contrast with display changes in the EKG, which generally 
15 requires a skilled cardiologist to evaluate. 

There are several possible factors which might cause 
a patient to transition from normal to VF. These factors 
may include drug overdose (especially overdose with an 
anti-arrhythmic which has a pro-arrhythmic effect in 
20 overdosage, e.g., quinidine . intoxication), excessive 
electrical stimulation, hypothermia, ischemia, and stress. 
In a preferred embodiment, a patient monitor may examine 
the patient's PPP so as to determine if the patient is in 
transition from normal to VF; this could indicate that one 
25 of these pro-arrhythmic factors is excessively present. 

The processor 103 may further process the PPP so as 
to measure the PPP's degree of deterministic chaos. 
Several techniques may be applied for this purpose: 

(1) The processor 103 may measure a Lyapunov 
30 exponent of the PPP. The Lyapunov exponent of the PPP is 
a measure of the degree to which nearby paths of the PPP 
diverge. The Lyapunov exponent is well-known in chaos 
theory and may be measured with available software. See, 
e.g., Wolf et al., "Determining Lyapunov exponents from a 
35 time series", Physica D 1985;16:285-317. 
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(2) The processor 103 may measure a fractal 
dimension of the PPP. The fractal dimension of the PPP is 
a measure of the degree to which the PPP forms a "space- 
filling" curve. The fractal dimension is well-known in 

5 chaos theory and may be measured with several techniques 
(e.g. correlation dimension or box-counting methods), for 
example as shown below: 

To measure the fractal dimension of the PPP, the 
processor 103 superimposes a rectilinear grid (comprising 

10 a set of boxes) on the PPP and counts the number of boxes 
which are cut by the PPP's trace. The processor 103 
varies the size of the grid and records each grid size and 
each count. The processor 103 then computes the constant 
k in the following relation: 

15 In (# of boxes cut) = k * In (# of boxes in gri(Q04) 

The constant k is a measure of the fractal dimension 
of the PPP. A value of k between about 3 and about 7, 
especially with a fractional component, implies that the 
PPP is likely to represent a process based on 

20 deterministic chaos, and therefore a patient who is close 
to (or actually in) VF. 

(3) The processor 103 may determine a Poincare 
section of the PPP and examine that Poincare section for 
indicators of deterministic chaos, as described herein. 

25 The processed PPP and Poincare sections may also be 
displayed for review by a human operator, whereupon any 
visible structure will be readily recognized. 

Figure 4 shows an example PPP 401 and a corresponding 
Poincare section 402. A Poincare section may comprise a 

30 line segment drawn across a part of the PPP. In general, 
such a line segment will be close to perpendicular to the 
trajectories of the PPP in a region of interest. 

The processor 103 may acquire the data points in each 
Poincare section or PPP and compute a statistical measure 

35 of anisotropy or inhomogeneity of those data points. One 
such measure is based on the mean and standard deviation 
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of those data points (these may be computed by statistical 
methods which are well-known in the art) . The ratio 

r - (standard deviation) / (expected value) (403) 
is a measure of the degree of clumping in the Poincare 
5 section. 

A greater value for r implies that the PPP is more 
likely to represent a process based on deterministic 
chaos, and therefore a patient who is close to (or 
actually in) VF. The value for r may be displayed for 
10 review by a human operator in comparison with a value for 
r for a normal patient, together with a set of confidence 
bands, as is well-known in the art, for indicating a 
degree of variation from a normal patient. 

The processor 103 may also compute other statistical 
15 measures of the Poincare section. 

The processor 103 may also determine a "time-lapse" 
Poincare section of the PPP. 

Figure 5 shows an example PPP 501 and a corresponding 
time-lapse Poincare section 502. A time-lapse Poincare 
20 section may comprise a set of data points selected from 
the PPP by selecting one data point every t seconds. The 
time-lapse Poincare section may be analyzed in like manner 
as the other Poincare section disclosed herein. 

A second aspect of the invention relates to detection 
25 and evaluation of heart disorders based on a frequency- 
domain transform of a patient EKG. 

Figure 6 shows a set of corresponding frequency- 
domain transforms, obtained by performing an FFT on the 
EKG signal. A first transform 601 corresponds to a first 
30 EKG signal (not shown) . A second transform 602 
corresponds to a second EKG signal (not shown) . 

In the first transform 601, representing a normal 
patient, the frequency spectrum shows that the energy of 
the corresponding EKG signal occurs primarily at a 
35 discrete set of frequencies. In the second transform 602, 
representing a patient exhibiting VF, the frequency 
spectrum shows that the energy of the corresponding EKG 
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signal has a continuous spectrum of frequencies, and has 
an energy peak 603. 

Part of this aspect of the invention is the use of 
both visual and mathematical techniques for analyzing 

5 frequency domain transforms, including for example 
calculation of a harmonic magnitude ratio (HMR) . To 
determine the HMR, a major peak or a central region of 
energy distribution in a spectrum of a frequency domain 
transform (such as an FFT) may be identified, and the HMR 

10 calculated as follows: A magnitude of the transform in 
the region of the identified point is determined (e.g., by 
summing the magnitude of the transform at the identified 
point and at surrounding points) , and is summed with the 
corresponding magnitude in the region of harmonic values 

15 of the frequency for the identified point. This sum is 
divided by a total magnitude of the transform for the 
entire signal; the ratio is defined as the HMR. 

. One method which is known for bringing a patient out 
of VF ("defibrillating") is to administer an electric 

20 shock across the patient's heart. This electric shock 
must generally have a substantial energy, e.g. 10-20 
joules, and may often cause tissue damage to the patient 
even if it is successful in defibrillating the patient. 
Multiple shocks may be required, generally of increasing 

25 energy. Accordingly, it would be advantageous to use a 
larger shock only when necessary, and it would be 
advantageous to use as few shocks as possible. 

Part of this aspect of the invention is the discovery 
that when the energy peak 603 of the frequency-domain 

30 transform 602 is at a relatively low frequency, a 
relatively low energy shock will generally suffice to 
defibrillate the patient. When the energy peak 603 of the 
frequency-domain transform 602 is at a relatively high 
frequency (also, when a secondary energy peak 604 appears 

35 in the frequency-domain transform 602 at a relatively high 
frequency) , it will require a relatively high energy shock 
to defibrillate the patient, if it is possible to 
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def ibrillate the patient by means of an electric shock at 
all. 

One application of this discovery is in automated 
implanted cardiac defibrillators (AICDs) , which attempt to 
5 automatically detect VP and to automatically administer a 
shock to def ibrillate the patient. 

Figure 7 shows an improved AICD 701. A patient 702 
is coupled to an AICD EKG 703, which acquires EKG signals 
and transmits them to an AICD processor 704, which 

10 controls a shock device 705 for administering a 
defibrillating shock to the patient 702. 

The improved AICD 701 also comprises (e.g., as part 
of the AICD processor 704) software for determining an FFT 
of the EKG signal and for determining the energy peak in 

15 that FFT. If the energy peak in that FFT is relatively 
low, the AICD processor 704 controls the shock device 705 
to administer a relativiely small shock to the patient. If 
the energy peak in that FFT is relatively high, the AICD 
processor 704 controls the shock device 705 to administer 

20 a relatively large shock to the patient, and may also 
signal an alarm 706 or other indicator that defibrillation 
may not be successful. 

A third aspect of the invention relates to detection 
and evaluation of drug toxicity based on a parameter time 

25 constant for an action-potential duration (APD) 
restitution curve or an action-potential amplitude (APA) 
curve which is constructed for the patient. 

Figure 8 shows a signal response of an individual 
heart muscle cell to a stimulus. This individual cell 

30 response is known in the art as "action potential". 

It is well-known in the art that a time duration for 
recovery 801 of an individual cell depends on factors 
including a resting period 802 which the cell has had 
prior to stimulus. It is also well-known in the art that 

35 an APD restitution curve can be constructed for a human 
patient with the use of an intracardiac catheter. 
However, the complete relation between the actual time 
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duration for recovery 801 based on the resting period 802 
is not known. 

Part of this aspect of the invention is the discovery 
that when the time duration for recovery 801 is plotted 
5 against the resting period 802 (diastolic interval) , the 
curve follows an exponential relation: 

APD = APD pl - A * exp(-DI/tau) (803) 
where APD pl is the plateau APD, A is a proportionality 
constant, DI is the diastolic interval, and tau is 
10 the parameter time constant 

The nonlinear nature of the APD restitution curve may 
promote deterministic chaos in response to excessive 
stimulus of the heart muscle cells. When the APD 
restitution curve is steeper (i.e., the parameter time 
15 constant tau is larger) , there is accordingly a greater 
predilection for the heart to enter VF. Thus, another 
peart of this aspect of the invention is the discovery that 
a normal patient will have a relatively low APD 
restitution parameter time constant, while a patient who 
20 is exhibiting drug toxicity (e.g., quinidine intoxication) 
will have a relatively high APD restitution parameter time 
constant. The restitution parameter time constant may 
also be used in monitoring cardiac stability, and in 
evaluating efficacy of anti-arrhythmic drugs. 
25 Experimental verification of the present invention 

has been achieved by the inventors. 

Experiment I. 

A mathematical study used PPPs, return maps, Poincare 
sections, correlation dimension, and spectral analysis to 

30 distinguish periodic, chaotic and random signals. PPPs 
were useful in distinguishing among all three classes of 
signals. Periodic signals showed clear, widely separated 
trajectories; chaotic signals showed banding, forbidden 
zones and sensitive dependence on initial conditions; 

35 random signals showed no clear internal structure. With 
the exception of noise effects, the only major difference 
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between the PPPs and the appropriately lagged return map 
was a 45 degree rotation. Poincare sections were also 
able to distinguish among the three classes of signals: 
periodic signals showed isolated points; chaotic signals 
5 showed ordered areas of apparent self-similarity; random 
signals showed a Gaussian distribution of points. 
Correlation dimension was more able to distinguish between 
chaotic and random signals than between chaotic and 
periodic signals. Spectral analysis using FFTs and 

10 harmonic magnitude ratio (HMR) was able to distinguish 
periodic signals , but were unable to distinguish between 
random and chaotic signals: HMRs of periodic signals were 
greater than 97%; HMRs of chaotic signals varied between 
17 and 80%; HMRs of random signals were approximately 40%. 

15 PPPs were greatly affected by noise, return maps were less 
affected, while spectral analysis was relatively immune to 
noise. It was concluded that PPPs, return maps, Poincare 
sections, correlation dimension and spectral analysis are 
all useful determinatives of chaotic systems. 

20 Experiment II. 

A mathematical study concentrated specifically on 
ability of spectral analysis to distinguish chaotic from 
random signals. In this experiment, two series of random 
signals were generated. The first series comprised 5000 

25 pseudo-random numbers which were smoothed using a method 
of least-squares approximation. The second series 
comprised white noise obtained from an analog-to-digital 
conversion board. Spectral analysis was performed by 
applying an FFT to the data, and searching for a broad 

30 band spectrum or a change from a narrow band to a broad 
band, which was presumed to be diagnostic of chaos. It 
was concluded that spectral analysis by itself was 
insufficient to unequivocally distinguish chaotic signals 
from random signals, and that additional tests such as 

35 PPPs and return maps were necessary for this purpose. 



WO 91/19452 



PCT/US91/04000 



14 

Experiment III. 

An experiment examined spectral analysis, 
visualization of PPPs and correlation dimension analysis, 
for usefulness in distinguishing between normal sinus 
5 rhythm and VF in dogs. Ischemia and re-perfusion were 
used as stress factors in closed-chest anesthetized dogs. 
Spectral analysis of the dogs having normal sinus rhythm 
revealed narrow-band spectra with fundamental frequencies 
at the sinus rate and harmonics extending beyond 50 Hz. 

10 PPPs were consistent with periodic dynamics, and dimension 
analysis revealed low dimensional behavior (1 — 2.5). In 
contrast, spectral analysis of the dogs having VF, 
revealed broad-band behavior with most of the energy at 6 
Hz, and with energy at all frequencies between 1 and 25 

15 Hz. PPPs showed constrained aperiodic behavior, and the 
dimensional analysis revealed higher dimensions (4-6) than 
that observed for the normal sinus rhythm dogs. Thus, all 
three techniques proved useful in distinguishing normal 
sinus rhythm from VF. 

20 Experiment IV. 

An experiment examined spectral analysis, visualiza- 
tion of PPPs, visualization of return maps, and 
correlation dimension analysis, for their usefulness in 
identifying VF in humans. These analytical techniques 

25 were applied to data from eight hypothermic patients 
undergoing spontaneous VF, and also to data from three 
normothermic patients with VF induced during electro- 
physiology testing. All patients had a broad band 
frequency spectrum (0-12 Hz), a low dimension (range 2- 

30 5) , and banding and forbidden zones on PPPs and return 
maps. It was concluded that spectral analysis, 
visualization of PPPs, visualization of return maps, and 
correlation dimension analysis are useful in detecting and 
evaluating VF. 
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Experiment V. 

An experiment examined spectral analysis, 
visualization of PPPs and correlation dimension analysis 
for their usefulness in distinguishing between normal 
5 sinus rhythm and VF in humans. VF in eight hypothermic 
human patients undergoing open-heart surgery was studied. 
In all patients, first and second order PPPs showed 
forbidden zones and banding, and an FFT revealed a 
relatively continuous power spectrum at all frequencies 
10 from zero to 25 Hz, with a majority of the power below 12 
Hz. In contrast , correlation dimension in all cases was 
less than 4 . It was concluded that multiphasic analysis 
of the data is preferable to reliance on a single 
analytical technique such as correlation dimension. 

15 Experiment yj. 

An experiment utilized spectral analysis and 
visualization of PPPs to elucidate the heterogenous nature 
of atrial fibrillation. In the experiment, the 
researchers induced acute fibrillation by a rapid train of 

20 stimuli to the atria of seven closed-chested dogs. PPPs 
based on the EKG data often inscribed well defined 
structures, and an FFT of the digitized EKGs showed peaks 
mostly below 15 Hz that were either discrete with clear 
harmonic components, or had continuous spectra that 

25 changed in a time- and site-dependent manner. It was 
concluded that both spectral analysis and visualization of 
PPPs are useful techniques for analyzing atrial as well as - 
ventricular fibrillation. 

Experiment V^. 

30 In an experiment, visual analysis of PPPs and the 

slope of an APD restitution curve were found to be useful 
for detecting and evaluating quinidine-induced VF in in 
vivo hearts. Quinidine was administered at 30 minute 
intervals over five hours, until either a total of 90-100 

35 mg/kg was administered or until ventricular tachycardia or 
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VF occurred, whichever came first. PPPs of the quinidine 
intoxicated cells demonstrated sensitive dependence on 
initial conditions and the presence of forbidden zones, 
and the corresponding FFTs showed continuous spectra. In 
5 contrast, PPPs of cells in a control dog were uniform and 
densely packed, and the corresponding FFTs showed discrete 
spectra. The initial slope of the APD restitution curve 
of quinidine intoxicated cells was much steeper, by at 
least an order of magnitude, than the slope of normal 
10 cells. It was concluded that quinidine toxicity 
correlates with the slope of the APD restitution curves. 

Experiment VIII. 

An experiment compared the slope of the APD and APA 
restitution curves with quinidine intoxication. Quinidine 

15 was administered (90-100 mg/kg) to eight dogs over a five 
hour period. Three untreated dogs served as controls. 
Ventricular and Purkinje cells from both treated and 
untreated dogs were then subjected to electrical 
stimulation with cycles from. 900 to below 600 msec. 

20 Shortening of the cycle length to 600 msec resulted in 
irregular dynamics of both APD and APA, including 
electrical alternans and bifurcation. The slope of an APD 
restitution curve was calculated, and found to be steeper 
in quinidine- intoxicated cells for both Purkinje fibers 

25 and ventricular muscle cells than the slope during 
quinidine washout or in normal untreated cells. The curve 
could be fit by the exponential equation given herein. 
APA changes were almost always correlated with the APD 
changes. In the three normal tissue preparations neither 

30 ventricular muscle cells nor Purkinje cells showed 
bifurcative behavior with respect to APD or AA. It was 
concluded that quinidine toxicity, and presumably other 
drug-induced pro-arrhythmic effects, correlate with the 
slope of both APD and APA restitution curves. 
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Experiment IX. 

In an experiment, quinidine- induced ventricular 
tachycardia and VF in dogs was analyzed using PPPs 
generated from action potential duration (APD) and action 
5 potential amplitude (APA) data. Both PPPs showed 
forbidden zones and sensitive dependence on initial 
conditions which are indicative of chaos. It was 
concluded that PPPs based on either APD or APA are useful 
in detecting and evaluating quinidine toxicity. 

10 Experiment X. 

In an experiment , EKGs of quinidine intoxicated dogs 
were analyzed by frequency spectra, phase plane plots, 
Poincare sections, return maps and Lyapunov exponents. In 
the control state and at therapeutic doses, PPPs were 

15 uniformly thick and showed no gaps, indicating that cycle- 
to-cycle variation was due to normal biological "noise". 
But as the quinidine dose was increased to intermediate 
levels (40-50 mg/kg) , PPPs showed clear non-uniform 
thickening, indicating sensitive dependence on initial 

20 conditions, and also showed marked banding (densely filled 
regions separated by divisions or gaps) . At these 
intermediate doses, Lyapunov exponents became positive and 
Poincare return maps also indicated nonrandom chaos. At 
still higher doses, PPPs became more complex. In two dogs 

25 that did exhibit VF (and not in another) there was a 
significant change in the PPP at the last pre-f ibrillatory 
dose: the development of a "funnel", a classic mechanism 
of chaos. Frequency spectra at all pre-f ibrillatory doses 
were discrete, with peaks at a fundamental frequency and 

30 multiple harmonics. It was concluded that chaos does 
occur during progressive quinidine intoxication, and that 
PPPs, and graphic and numeric analysis based on the PPPs, 
are better indicators of chaos than frequency spectra. 



WO 91/19452 



PCT/US91/04000 



18 

Experiment XI. 

In an experiment, quinidine toxicity in dogs was 
analyzed using PPPs generated from APA and APD data. EKG 
recordings were made at various driving rates from 1000 to 
5 500 msec. Increase in the driving rate from 1000 to 500 
msec caused the progressive appearance of higher order 
periodicities (period 3 and 4), Phase locking was seen 
with a stimulus (S) response (R) pattern repeating 
periodically in all 4 preparations at S:R ratios of 2:1, 

10 5:3 f 3:2. At faster drive rates aperiodic variations in 
APA and APD were observed. A number of intermediate 
stages that presage chaos were also seen in the quinidine 
intoxicated fibers. These results further demonstrate the 
usefulness of the methods of the present invention to 

15 detect both quinidine intoxication and precursor stages to 
intoxication. 

Experiment XII. 

In an experiment, quinidine toxicity in dogs was 
analyzed using PPPs generated, from APA and APD data. 

20 Electrical stimuli were used to drive cardiac tissue at 
various rates from 2000 to under 300 msec. These stimuli 
caused steady alternans (bifurcation) in APD and APA of 
108 ± 36 msec and 12 ± 9 millivolts respectively. Further 
increase in driving rates gave rise to irregular dynamics. 

25 This transition was preceded by various repeating 
stimulus-response ratios (phase-locking) for up to fifty 
consecutive beats. No such dynamics could be induced in 
three non treated (control) tissues. The APD restitution 
curve had significantly (p < 0.05) steeper slope than six 

30 control fibers. stimulus-response latency remained 
constant at 6-9 msec. PPPs of the APDs during the 
irregular dynamics showed sensitive dependence on initial 
conditions and forbidden zones consistent with chaos 
theory. These results further demonstrate the usefulness 

35 of the methods of the present invention to detect both 
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quinidine intoxication and precursor stages to 
intoxication. 

Experiment XIII, 

An experiment used spectral analysis , PPPs, Poincare 
5 sections, Lyapunov Exponents and dimension analysis to 
analyze computer simulated waveforms including sine waves, 
modulated sine waves, square waves, saw toothed waves, and 
triangular waves. The researchers added random noise to 
the waveforms at 1%, 10% and 20%. The experiment further 

10 used the same analytical techniques on EKG data from 
anesthetized dogs in which VF was precipitated by five 
different interventions: quinidine intoxication; premature 
electrical stimulation followed by quinidine intoxication; 
coronary occlusion; reperfusion of acutely ischemic 

15 myocardium; and global hypothermia. The preliminary 
results showed that PPPs and Poincare sections in dogs 
undergoing ventricular fibrillation were consistent with 
chaos, while spectral analysis was not suggestive of 
chaos. The researchers concluded in part that VF can be 

20 described as chaotic electrophysiological behavior, but 
that single methods of analysis are not sufficient to 
detect such behavior. 

One conclusion which may be drawn from the research 
cited herein is that the analytical value of each of the 

25 aspects of the invention may be enhanced through 
combination with one or more of the other aspects of the 
invention. A preferred embodiment of the present 
invention may include a combination of the aspects of the 
invention described herein. One preferred embodiment may 

30 comprise multiphasic analysis of a PPP (e.g., visually 
with a display, graphically with Poincare sections, and 
numerically with Lyapunov exponents and correlation 
dimension) , frequency spectral analysis, and mathematical 
analysis of an APD restitution curve. 
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Alternative Embodiments 

While preferred embodiments are disclosed herein, 
many variations are possible which remain within the 
concept and scope of the invention, and these variations 
5 would become clear to one of ordinary skill in the art 
after perusal of the specification, drawings and claims 
herein. 

It would also become clear to one of ordinary skill 
in the art that embodiments of the invention may comprise 

10 means for continuous monitoring of drug toxicity, atrial 
fibrillation, ischemia or other heart conditions, such as 
during surgery or patient recovery from surgery. 
Moreover, embodiments of the invention may comprise means 
for indicating heart conditions which sure detected to 

15 attending medical personnel or to the patient. In one 
preferred embodiment of the invention, means may be 
provided for directing the patient (when a heart disorder 
is detected) to contact a physician or to proceed to a 
nearby hospital for treatment. 
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Claims 

1. A method for detecting heart disorders, 
comprising the steps of 

receiving an electrocardiogram signal; 
5 determining a phase-plane plot of said 

electrocardiogram signal; and 

determining if said phase-plane plot indicates 
a heart disorder. 

2. A method as in claim 1, wherein said heart 
, 10 disorders are one of the group: ischemia, electrical 

instability, drug toxicity. 

3. A method as in claim 1, wherein said phase- 
plane plot comprises a multi-dimensional plot of at least 
two variables. 

15 4. A method as in claim 1, wherein said phase- 

plane plot comprises a multi-dimensional plot of at least 
three variables. 

5. A method as in claim 1, wherein said phase- 
plane plot comprises a plot of a signal voltage against a 

20 derivative of said signal voltage. 

6. A method as in claim 5, wherein said derivative 
is a first derivative. 

7. A method as in claim 1, wherein said step of 
determining if said phase-plane plot indicates a heart 

25 disorder comprises the step of determining a fractal 
dimension of said phase-plane plot. 

8. A method as in claim 1, wherein said step of 
determining if said phase-plane plot indicates a . heart 
disorder comprises the step of determining a Lyapunov 

30 exponent of said phase-plane plot. 
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9. A method as in claim 1, wherein said step of 
determining if said phase-plane plot indicates a heart 
disorder comprises the step of determining the presence of 
at least one of: banding, forbidden zones, nonuniform 

5 thickening, periodicity, aperiodicity. 

10. A method as in claim 1, wherein said step of 
determining if said phase-plane plot indicates a heart 
disorder comprises the step of determining the presence of 
multiple alternants. 

10 11. A method for detecting heart disorders, 

comprising the steps of 

receiving an electrocardiogram signal; 
determining a phase-plane plot of said 
electrocardiogram signal; 
15 determining a Poincare section of said phase- 

plane plot; and 

determining if said Poincare section indicates 
a heart disorder. 

12. A method as in claim 11, wherein said phase- 
20 plane plot comprises a multi-dimensional plot of at least 

two variables. 

13. A method as in claim 11, wherein said phase- 
plane plot comprises a multi-dimensional plot of at least 
three variables. 

25 14. A method as in claim 11, wherein said phase- 

plane plot comprises a plot of a signal voltage against a 
derivative of said signal voltage. 

15. A method as in claim 14, wherein said derivative 
is a first derivative. 
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16. A method as in claim 11, wherein said step of 
determining if said Poincare section indicates a heart 
disorder comprises the step of determining a statistical 
measure of at least one of: anisotropy, inhomogeneity. 

5 17. A method for detecting heart disorders, 

comprising the steps of 

receiving an electrocardiogram signal; 
determining a frequency domain transform of said 
electrocardiogram signal; 
10 determining if said frequency domain transform 

indicates a heart disorder. 

18. A method as in claim 17, wherein said step of 
determining a frequency domain transform comprises the 
step of performing a Fourier transform. 

15 19. A method as in claim 17, wherein said step of 

determining if said frequency domain transform indicates 
a heart disorder comprises the step of determining if said 
frequency domain transform comprises a continuous 
spectrum. 

20 20. A method as in claim 17, wherein said step of 

determining if said frequency domain transform indicates 
a heart disorder comprises the step of determining if said 
frequency domain transform comprises an energy peak at a 
relatively high frequency. 

25 21. A method as in claim 20, wherein said relatively 

high frequency is greater than about 6 Hz. 

22. A method as in claim 20, wherein said relatively 
high frequency is greater than about 10 Hz. 

23. An automated implantable cardiac defibrillator, 
30 comprising 
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means for receiving an electrocardiogram signal 
from a patient; 

means for administering a def ibrillating shock 
from said patient; 
5 means for processing said electrocardiogram 

signal and for controlling said means for administering, 
wherein said means for processing causes said means for 
administering to administer a shock whose energy depends 
on a frequency domain transform of said electrocardiogram 
10 signal. 

24. A defibrillator as in claim 23, wherein said 
energy is relatively greater when said frequency domain 
transform has an energy peak at a relatively greater 
frequency. 

15 25. A method for detecting drug toxicity for a 

patient, comprising the steps of 

constructing a relation between a diastolic 
interval and an action potential duration for said 
patient; 

20 determining a value of a time constant for said 

relation; and 

determining if said value indicates drug 

toxicity. 

26. A method as in claim 25, wherein said step of 
25 determining if said value indicates drug toxicity 

comprises the step of determining if said value is 
substantially above a normal value. 

27. A method for evaluating effectiveness of a drug 
for a patient, comprising the steps of 

30 constructing a relation between a diastolic 

interval and an action potential duration for. said 
patient; 
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determining a value of a time constant for said 
relation; and 

determining if said value indicates that 
administration of said drug to said patient has an anti- 
5 arrhythmic effect. 

28. A method for evaluating effectiveness of a drug 
for a patient, comprising the steps of 

determining a set of action potential duration 
restitution data for said patient; 
10 constructing a phase plane plot from said data; 

and 

performing a multiphasic analysis of said phase 
plane plot to determine if said phase plane plot indicates 
that administration of said drug to said patient has an 
15 anti-arrhythmic effect. 

29. A method as in claim 28 , wherein said 
multiphasic analysis comprises at least one step of 
displaying said phase plane plot, computing a Poincare 
section of said phase plane plot, computing a Lyapunov 

20 exponent of said phase plane plot, or computing a 
correlation dimension of said phase plane plot. 

30. A method for detecting heart disorders, 
comprising the steps of 

receiving an electrocardiogram signal; 

25 performing a plurality of the following steps to 

determining if said electrocardiogram signal indicates a 
heart disorder: (a) performing a multiphasic analysis of 
a phase-plane plot of said electrocardiogram signal, (b) 
performing a spectral analysis of said electrocardiogram 

30 signal, or (c) performing an analysis of an APD 
restitution curve computed in response to said 
electrocardiogram signal. 
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31. A method as in claim 30, wherein said step of 
performing a spectral analysis comprises the step of 
calculation of a harmonic magnitude ratio. 



WO 91/19452 



PCT/US91/04000 



1/4 





/ 



70/ 



ACID 




alaam\ 







706 




SUBSTITUTE SHEET 



WO 91/19452 



PCT/US91/04000 




SUBSTITUTE SHEET 



WO 91/19452 PCT/US91/04000 

3/4 




SUBSTITUTE SHEET 



WO 91/19452 



PCT/US91/04000 




